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Abstract 

The issue discussed is a thermodynamic version of the Bern-Kosower master amplitude 
formula, which contains all necessary one-loop Feynman diagrams. It is demonstrated how 
the master amplitude at finite values of temperature and chemical potential can be formu- 
lated within the framework of the world-line formalism. In particular we present an elegant 
method how to introduce a chemical potential for a loop in the master formula. Various 
useful integral formulae for the master amplitude are then obtained. The non-analytic prop- 
erty of the master formula is also derived in the zero temperature limit with the value of 
chemical potential kept finite. 
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1. Introduction 



Based on the field theory limit of string theory (with infinite string tension limit), a very 
elegant method was invented several years ago: the Bern-Kosower (BK) rules to obtain one- 
loop gluon scattering amplitudes in a compact form |jj (see also 0). For example, five gluon 
scatterings are efficiently calculated by using these rules []|] , and various field theory limits have 
been studied along the line of the BK formalism: perturbative gravity (4j, super Yang-Mills 
theory bosonic string theory approach ||, and multi-loop generalizations ||, ||]. 

The most important and conspicuous point in the BK formalism is that all Feynman dia- 
grams are included in a single master integrand. In this formalism, we hence do not compute 
the loop integrals and the Dirac traces of respective Feynman diagrams. These facts can also 
be achieved by another approach, called the world-line formalism, which reformulates Feynman 
(field theory) amplitudes similar to string theory amplitudes; i.e. the field theory amplitudes 
can generally be obtained as a path integral average of vertex operators [pX|1-p7|. 

In fact, there have been established many examples: (j) 3 theory g(| ITT), QED gj, |l|, axial 
vector and pseudo-scalar couplings [14, 15], Yang-Mills theory |L6|, [17| (see also |l^, [OJ), and 



more references can be found in [20]. Both world-line and string theory methods help each 
other, and are useful to get an insight for solving mutually related problems in each own way; 
for example, this viewpoint has been very useful for specifying pinching limits and corners of 
string moduli in the multi-loop analysis ]|]. It is very interesting that these two methods, which 
entirely differ from the conventional Feynman diagram technique, improve the computational 
efficiency for obtaining Feynman amplitudes. 

However, compared to these developments, their thermodynamic versions have not been 



studied very much [22|-|2(| from the viewpoint of general formulation for constructing a master 
amplitude formula. In particular, there has been no general and convenient method how to 
introduce a chemical potential for a loop in the world-line formalism. Before stepping in an un- 
explored calculation by using the formalism, it is important to establish a definite and universal 
foundation in the first place. In this paper, therefore as a basic step on the thermodynamic 
generalization of world-line field theory, we present a universal and fundamental prescription 
of one-loop iV-point amplitudes at finite values of temperature and chemical potential without 
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help of any standard calculation. We shall study the amplitude of a particular form (the master 
formula), where the loop integration and the Matsubara summation are a priori finished and 
Feynman's parameter integrals are only left — laying emphasis on the point that we never 
mean the Feynman integrals as a single Feynman diagram, but as a sum of all diagrams. It is 
certainly nontrivial to introduce a chemical potential (as well as a temperature) with keeping 
this advantageous point pertaining to the master formula intact. 

We address the following points in the matter of the thermodynamic generalization (at finite 
values of the temperature and the chemical potential [i for a loop). First, we present a 
formulation of the thermodynamic amplitudes along the same line as the non-thermodynamic 
world-line formalism. In particular, the way of introducing the chemical potential is a non- 
trivial problem. Since we do not introduce any idea of the continuous or discrete momentum 
integration/summation, we have to find an alternative to the shift of internal discrete momenta: 

<jj n — >■ UJ n + in . (1.1) 

This situation may be understood in the following way: In the standard method, the inverse 
temperature j3 is introduced by summing up all topological different S 1 paths along the zeroth 
component (imaginary time) direction. On the other hand in the world-line formalism, this 
procedure modifies the path integral of a corresponding periodic world-line field x (<r); < a < 
1, into a summation of the path integrals with x°(a) shifted by 

x°{a) -► x°{a) + nfia . (1.2) 

Although the radius (3 of S 1 and the world-line circumference (the unity) have nothing to do 
with each other, the shift (1.2) involves the world-line coordinate a as well. In this sense, the 



way of prescribing the temperature and hence of the chemical potential become different from 
the standard Matsubara formalism. To introduce the chemical potential, one may of course 
transform a world-line amplitude formula to the Feynman-Matsubara form, and then resume 
the world-line form after some efforts to apply the shift ([O]). However, such a calculation 
does not utilize the merits of the master formula at all, and there should be a more direct and 
transparent method to introduce fx within the world-line formalism (without tracing back and 
referring any internal loop calculation). 



To this end, we propose a new rule, instead of to introduce the internal chemical 

potential. It is simply done by applying the new shift procedure 



u — ► uj + i/j, , (1-3) 

where Co is an average of the zeroth components of external continuous/discrete momenta k®\ 
j = 1, 2, • • • N, with the summation weights aj (the local coordinates of the external legs on 
the closed loop). The prescription ( |1.3[ ) is neither conceivable nor explicable from the standard 
method (|1.1| ), because uo n is the internal momentum while Co concerns the external one. The 
parameter uj will easily be identified in due course, if we adopt a statistical parameter to discern 
between boson and fermion loops when summing up the S* 1 paths: This parameter dependence 
should vanish at the [3 — ► oo limit as expected in the non-thermodynamic world-line formalism. 

Another non-triviality in this formalism is how to derive a non-analytic property at (3 = oo 
with finite \x from the master amplitude formula. Since the non-analytic property can be 
derived from the (3 = oo limits of pure thermodynamic parts, we first have to separate a pure 
thermodynamic part f jy from a full thermodynamic amplitude T^ 1 . If the master integrand 
of r^y 1 is composed only of the Jacobi 6-function, the story is simple as expected. However, 
in a more complicated case like a photon scattering, the master integrand is not such a simple 
form but a product of a (3- and //-dependent operator Vp^ and the Jacobi G-function part Kp^- 
The pure thermodynamic parts of these quantities, and ICp^, can easily be separated from 
their original full quantities, but we show that the pure thermodynamic part f is non-trivially 
given by Vp^ x fCp^ against naive expectation. (The simple arithmetical splittings indicate one 
more contribution, however we shall show that it vanishes). Separating the f ^ in this way, we 
then analyze the non-analytic property of the master amplitude at zero temperature. 

This paper is organized as follows. In Section we explain our notations, definitions and 
the general structure of master amplitudes at zero temperature. In Section |3], we derive a set 
of general formulae for a master amplitude at finite (3 and /U parts by parts: path integral 
normalization and (scalar) kinematical factor in Section [O], and reduced kinematical factor in 



Section 3.2. In Section |3.3|, we show a general master formula for the purely thermodynamic 



part of the full amplitude T^i an d also prove the above statement, i.e., f jy ~ Vp^ x Kp^. 



The replacement rule (|l.3| ) is verified in Appendices A-C from the viewpoints of both Feynman 
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rule's calculation and the world-line formalism. In Sections |] and ||, for arbitrary N, we derive 
various integral formulae by analyzing two kinds of O-function representations, and check their 
consistency. Several explicit results are also presented up to N = 5. In Section ^, we derive the 
non-analytic property of the master formula in the zero temperature limit with the chemical 
potential kept finite. Section [7| contains summary and conclusions. 



2 Notations and definitions 

First, we summarize the general structure of the one-loop iV-point master amplitudes 1 in 
the non-thermodynamic world-line formalism. We refer the reader to Refs. [12, [16|, [2(| for more 
details. The master amplitude of general form (for N external momenta j = 1,2, ••• , N; 
(j, = 0, 1, • • • , D— 1) is written in terms of the closed path integrals of bosonic x^{a) and fermionic 
ip^{a) world-line fields as follows: 

T N = - T— <f > Vx»{o)V^(<j)e-ti C[a)d<T TT Vj (2.1) 
2 Jo s J 

with the world-line Lagrangian and the vertex operators 

V 3 = s I' davAxial^a)^^ ; j = 1, 2, - • • , N , (2.3) 



JO 

where kj ■ x stands for the Lorentz contraction, and we often omit the Lorentz indices as long 



as obvious. The zero mode integral of the bosonic path integral should be excluded |18| . The 
explicit form of Vj depends on what particle is inserted in the loop as an external leg; for 
example, Vj = 1 for 3 theory, and is in Eq.( p.2|) for the photon vertex case. Note that our 
world-line coordinate a is dimensionless, and is related to the standard notation |20| by scaling 
r = sa. For the path integral average of a general quantity F, 

< F(x, tp) >= AT' 1 j VxVi>e~ C{a)dr7 F{x, ip) , (2.4) 

one may use the Wick contractions with 

< x^jpi) x> 2 ) >= -sg^G(a 1 ,a 2 ) , (2.5) 

1 We assume that a particle change does not occur while circulating along the loop. 
2 The Euclidean metric is given by g^ v = -fr LV . 
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< naxWfa) >= ^signiai - a 2 ) , (2.6) 

where M is the path integral normalization 

M = j> VxV^e' So C{a)da = e~ sm ' 2 (4tts)- d/2 , (2.7) 

and G the bosonic world-line correlator 

G(ai,<jj) = |crj — (Tj\ — (cjj — cr,) 2 = Gij . (2.8) 

Performing the path integrals (or Wick contractions), we arrive at the following master 
amplitude formula: 

rl rl rl /*oo 

T N = c dux da 2 --- da N dss N ~ x VxK, (2.9) 
Jo Jo Jo Jo 

where the constant c depends on a theory |^ (particle circulating in a loop); For example, 

1 for a (neutral) scalar loop 

2 v ; p 

— ^tr[l] for a fermion loop, 

where the tr[l] expresses the trace of a unit matrix in the .D-dimensional gamma matrix space. 
For a gluon loop, c is no longer a simple constant (T(^, The quantity /C, which we shall call 
the kinematical factor (multiplied by the normalization jV), is defined by the world- line path 
integral average of N scalar vertex operators [10, 16 



N , N 



K = N < f[ e ik * x ^ >= Ivx^V^ (f[ e ik i- x ^)e~ So C{a)da . (2.11) 



The V stands for an effective vertex function, which is obtained by 



< njli v j exp[ikj ■ x((7j)] > 
< IljLi ex P[^i " x {°~j)\ > 



V= J , - / -V ■ (2.12) 



We shall refer to this quantity as the vertex structure function/ operator, which corresponds to 
the quantities called the reduced kinematical factor ]l]] or the generating kinematical factor [|l(| . 
At this stage, the V is still a function of s and o~j (j = 1, 2, • • • , N), however it will be generalized 
to an operator in the thermodynamic case. 

3 It is also related to an over-counting factor of a-integration regions j2lfl. 



If the vertex structure function can be expanded in the form 



V = Y / ais l ~ N exp[-sbi] , (2.13) 

where the coefficients ai and b\ may not depend on s, but on aj, the amplitude ( fO| ) is obtained 
as the sum of the 'partial' amplitudes 

Ai = / dss l - l IC = / dss^N < TT e ik r x ^) > , (2.14) 
Jo JO jj[ 

with shifting m 2 to m 2 + fy. 

In the case of finite /?, we assume the zeroth components of the external momenta k^ to be 
the bosonic Matsubara frequencies [] 

k °j = "kj = -jnj , (2-15) 

as well as 

kj ■ x = uJkjX — kj ■ x . (2-16) 
We use the notations for the counterparts of Eqs. (g]|), fl2~Tl| ), ( j^) , O), (|2.14|) 



T P N , Kp, N/3, Vp, A?. (2.17) 
In the case of finite fi and j3, we denote 

I$\ A//3 M , V^, ./if". (2.18) 

We explain how to define these quantities in the next section. 



3 The thermodynamic generalization 

3.1 The kinematical factor JCp^ and the normalization A//? M 

We take a two-step procedure to introduce the temperature and the chemical potential. As 

the first step, we consider the fi = case. Basically we follow the same method as discussed in 
4 We do not consider the fermionic external states, since the bosonic external states are only well-formulated 
in the world-line formalism. However one may formally generalize to the fermionic states. 
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Refs. |23[]-[26|, and we start with the following slightly general definition: 



CO 



= > e £ Tn , (3.1 



n=— oo 



where e is a constant related to the statistics of the loop. For example, the e = 2 case 
corresponds to a fermion loop, and the e = oo case to a bosonic loop. Without specifying the 
value of e, we deal with both cases simultaneously (formally fractional statistics as well). The 
Ttv on the right-hand side (r.h.s.) of the above formula denotes the path integral representation 
(2.1), and the replacement (|1.2|) should be applied to the x° fields for all variables aj. In a 
nutshell, we have only to replace the bosonic path integral in the following way: 



oo 

2mvi 



Vx^ -> & Vx^ = V e~ & Vx» . (3.2) 

For simplicity, let us consider the V = 1 case, or the 'partial' amplitude An (the I = N 
term in ( 2.13[ )). In this case, we have only to generalize the kinematical factor K, to the finite 
temperature version JCp; 

oo N 

def. y, e *ptf < TTjk j .x{* i ) > (3 3) 

n=— oo j=l 

= (4-)^e-^ 2 9 3 (i + §,^), (3.4) 
where the definition of ©3(1", r) is 

00 

Q 3 (v,t) = J2 e n2 ™e 2nvni , (3.5) 

n=— 00 

and we have introduced the following two s-independent quantities: 

N 

M 2 = m 2 - k i ■ k o G i3 > ( 3 - 6 ) 

i<j,=l 

N 

Q = ^2ajU kj . (3.7) 



Here the two-point correlator (world-line Green function) Gij , given by Eq. ( |2.8| ) , looks different 
from the one used in Ref . |p3|| , however the value of M 2 does not differ under the condition of 
momentum conservation regarding the external legs. The sign of M 2 seems not always to be 
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positive in spite that < < 1/2. We then assume M 2 > in the following discussion, 
choosing the off-shell symmetric point (satisfying the momentum conservation constraint) 

kikj = 5ijk 2 + (Sij - 1)^Y^ 2 • (3-8) 

For the V / 1 case, we have to apply Eq.([T^) to the V and tC parts at the same time, however 
the things are rather straightforward. We refer the reader to Appendix C for an example. 

The second step is to introduce the chemical potential for the internal loop. One may do it 
with applying the Jacobi transformation to the quantity KLp; i.e., rewriting Eq. fl3,4| ) in such a 
way to revive the discrete summation over internal Matsubara frequencies u n , one may perform 



the replacement (|1.1| ) (see [25] for more details). However, this is a roundabout way, since 
one has to make the Matsubara summation re-appear in spite of dealing with the world-line 
formulation, where the integration and summation of the loop are already performed. Instead, 
we propose a much simpler and direct alternative method to steer clear of this problem. It can 
be promptly done by the replacement fll.3j ), 

u> — ► ft = u) + i/jl . (3-9) 

This shift looks similar to the one (|l]l]) , however we stress that our shifting parameter is not the 
internal frequency but an average of external ones (cf. Eq.( |3.7D ). In this sense, this prescription 
is nontrivial. For a rigorous reader, we put a justification of this replacement in Appendix A 
(the case of V = 1), and also refer to Appendix B for more complicated case (the V / 1 case). 
After all, the general path integral average for finite [i can be obtained in a couple of simple 
replacements: 

Ivx" -» Ivx^ -> I = / Vx^ . (3.10) 

J J/3 Jf3n J/3 ^ n 

Note that the fermion path integral does not change by any means. 

Let us write down the thermodynamic kinematical factor K/3n and the normalization factor 
Npu- Applying the above replacement to Eq. fl3,4| ), we thereby yield the desired thermodynamic 
kinematical factor for finite (5 and ji: 

K Pli d ^ Eq.(3.4) _ (3.11) 

= (4^)^e-^ 2 G 3 ^,^), (3.12) 



where we have introduced the shorthand notation [i by reason of analogy to the vacuum am- 
plitude of a scalar loop 

27TZ 2-7T 

p, = fi-iui - — = -*(fi + — ) . (3.13) 



We refer to the representation ( 3.1 2| ) as the first representation. This representation with [i = 



(fj, = oj = 0, e = oo) is studied in Ref. [p3| . We also obtain another expression (which we shall 
call the second representation) through the Jacobi transformation, 

^ = I(4 OT )W'<«^>e3(^,^). (3.14) 

Note that for a fermion loop (e = 2), the above first and the second G3 representations become 
the 84 and 62 representations respectively. 

We can extract the thermodynamic quantity Afpu corresponding to the path integral nor- 
malization ( |2.7| ) from the kinematical factor JCp^. If we rewrite Eq.( |3.12| ) as 

= e- sm2 (47r S )- D / 2 e 3 (^,^)exp[ S f;A: J .^G 4J ] 

= ^4^)^e-^ 2 - 2 )e3(^,^)<ne^^)>, (3.15) 

the non-correlator part can be regarded as an overall normalization to the iV-point correlator of 
zero temperature type (the zeroth components k® are formally regarded as continuous variables 
here). Paraphrasing this fact in analogy to Eq.(2.11) 



N . j N 

= Afp„ < JJ e ik i< a ^ >= f VxViie-Jo c(a)da J[ e ik r<°i) ; ( 3 .i 6 ) 
3=1 J ^ 3=1 

the thermodynamic version of the path integral normalization (for finite (5 and y) is found to 



be 



^ = ^(4^ e -»^)e 3 (^,^). (3.17) 

This is exactly the same normalization as assumed in Ref. |22j] (taking a mass term inclusion 
into account). 



3.2 The vertex structure operator V, 



In this subsection, we discuss the part of vertex structure (V / 1). Because of diversities of 
explicit forms of V, there is no concrete formula such as Eq. (|3.12D , However, we can derive a 



general property of the pure thermodynamic part Vgu- To explain this, we classify V into two 
categories by the criterion whether or not V contains local world- line variables Gj. 

First, let us start with the first category, the u-independent V case. Apparently, V = 1 is 
the case. A nontrivial example of this category is the ir° — ► 27 decay (e = 2, D = 4) without 
background field [25| ]: 

cV = -tr[l]m\e 2 e flupa e l te%k^kZ , (3.18) 

where A is the pseudo-scalar coupling, m and e the (space-time) fermion mass and the QED 
coupling constant. Another nontrivial example is the effective potential of a fermion loop 
(e = 2, N = 0) in a constant magnetic field in D dimensions |p2| . In this case, V depends on 
the integration variable s, 

cV = —=br[l] sBcoth(sB) , (3.19) 



and this is the case of the expansion ( 2.13j ), if we use 



coth(sB) = 1 + 2J2 e~ 2nsB • (3-20) 



n=l 



With the shift M 2 — > M 2 + 21B, this case is essentially described by the 'partial' amplitude 



A\ defined by Eq.( 2.14 ). Now, at finite values of j3 and /i, we can, in principle, express any 



thermodynamic vertex structure function Vp^ as 

V/^ = V + V 0IX , (3.21) 

where Vp p denotes the purely thermodynamic part. However, in this category, as can be seen 
from the above examples, we simply have 

= V , V Pli = (for cr-independent V) , (3.22) 

because the V's of this category do not contain the local world-line variables aj, strictly speaking 
the Gij which are generated by the bosonic field correlation from Vj — Recall that the (5- 
dependence only appears through the shift x° — > x° + nf3a in the formula (|3.1|). The in- 
dependence can not be created from the quantity which does not contain a x field correlation. 

Second, let us consider the other category, the cr-dependent V case. In this category, we 
obtain a nonzero structure function Van,- We here consider the photon polarization case by way 
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of example. In this case, the V is given by Eq.( C.5 ) at zero temperature [16]: 

1 



V = ei • e 2 -Gi2 + ei • k 2 e 2 ■ h(Gi 2 ) 2 + e\ ■ e 2 k\ ■ k 2 — e\ ■ k 2 e 2 ■ k\ , (3.23) 



s 



where dots on G\ 2 means the first and the second derivatives w.r.t. the first argument of G± 2 . 



The thermodynamic generalization can be done by applying the formulae ( |3.1|) and (3J3). For 
ease of presentation, we put the computational details in Appendix C, and a comparison with 
the Feynman diagram technique is also in Appendix B. For finite (3, we derive Eq.( p.l^ ): 

Vfi = V - efc2(~) 2 " ■ h - ■ k 2 )G 12 , (3.24) 



and then shifting u> — > Q, we acquire the operator part given by Eq.(C.19): 



= - £ o e o(^) 2 " \^&* ■ h - eg ei • k 2 )G 12 . (3.25) 

Note that the origin of d/dVt is the Wick contractions of the bosonic world-line fields furnished 
in the Vj parts of photon vertex operators (see Eqs.(C2) and ( C.10| )), and it always happens if 



a vertex operator comprises a bosonic field in such a way. 

We therefore conclude that if V includes Gjj or its derivatives, Vp^ gives rise to a differential 
polynomial in d/d£l, and otherwise Vp^ = 0. An important result followed from this fact is 

V P(t x K = (for all V) . (3.26) 

This is a model-independent result, and we shall make use of this relation in order to decouple 
the pure thermodynamic part f^f from the total ./V-point amplitude T^ 1 . 

3.3 The master formulae fjf and A^ 

Gathering the formulae obtained in the above sections, we compose the purely thermody- 
namic master amplitude f jJf, and define the thermodynamic 'partial' amplitudes Ai ■ Applying 
the following formula to the first representation fl3.12|) 

oo 

e 3 (u, r) = 1 + 2 e n2 ™cos(2n™) , (3.27) 

n=l 

we separate the pure thermodynamic part /Cg^ from JCp^ as 

fcpn =1C + ICpn , (3.28) 
11 



where 



£^ = 2(4*s) — e 



2_, e is cosh(n/?//) , 

n=l 



or for the second representation ( 3.14 ) 

2 , l-D 



n=l 







For a given Vp^, the thermodynamic master amplitude T^ 1 is calculated as 



def. 



N 



c I dai / do2 

J 



da N I dss N 1 Vfl„ x /C/?„ . 
o j o 



Using the decompositions ( p.21| ) and (3.28) with the general formula ( 3.26j ), we find 



where 



\Pn def. 



N 



L Af ~~ 1 iV + 1 at j 



1 /"l /*1 /"OO 

c / do"! / dcr 2 • • • / dcrjv / dss N ~ 1 Vpn x /C,g M . 
j JO J 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



Therefore we have only to separate the pure thermodynamic part of fcpu in order to obtain 
the purely thermodynamic part • no separation in the Vp^ part at all (Note the difference 
between Eqs.( |3,31 ) and (3.33)). 

As can be inferred from the examples in Section [T^, a general expansion form of Vp^ is 

J-N ° 



V, 



Pfi 



a lr> s ' 



LneZ 



exp[-s6i] . 



(3.34) 



where the coefficients a\ n and b\ may not depend on s, but on Gj (j = 1, 2, • • • , N). In analogy 
to Eq,( [2.14D , the relevant 'partial' amplitudes for Eqs.( 3.31 ) and (|3.33| ) are defined as 



A 



Pfi def. 



dss 1 1 ICp fl 



m 2 ^m 2 +bi ' 



(3.35) 



7P/i def. 



ds J~ x K f 



m 2 — >TO 2 +6i 



(3.36) 



Ai - j us a r^p ( , 
and of course the following relation holds: 

= Ai + A^ . (3.37) 

In the following sections, we present various integral representations for the pure thermody- 
namic parts of the iV-point 'partial' amplitudes without specifying any values of N, D and 
e. These parts are the essential quantities to analyze the zero temperature limits with revealing 
the non-analyticity on ji. 
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4 The integral formulae from the first representation 

In this section, we derive various integral formulae for A.} based on the first representa- 
tion ( |3.12 ). This representation is examined in the special case I = 1 (with N = and e = 2), 



and actually (4-7r) D / 2 .A^ is the function Op analyzed in Ref.[22|: 

Op{m) =4^(-ircosh(n/? / u)(^) 1 - D / 2 ^ /2 _ 1 (n/?m) . (4.1) 

n=l 

We want to generalize this function to more generic I and e cases, and it is convenient to mimic 
the computational technique of Ref.(2^], introducing the parallel notation 

0f\M) d =- (4vr)f^ ; k = 2l + \-D. (4.2) 

Hereafter, for simplicity we set 

h = , (4.3) 

and for later convenience, we also define 

d = 3-k = D + 2-2l. (4.4) 
Prom Eqs.(^29|) and ( |336|) ,the pure thermodynamic part is now given by 

POO 2 2 

0f\M) = 2(4vr)T / ( f ss / - 1 (47rs)^e- sM2 £ cosh(nf3fi) . (4.5) 

•* n=l 

Performing the s-integration, we obtain 

oo n 

0f{M) = 4 £ co Mnm{^) l - d,2 K d/2 ^(n(3M) , (4.6) 

n=l 

where K u (z) is the modified Bessel function of second kind. This is a generalized version (for 
arbitrary e and /) of the above function Op. The aim is to obtain integral representations with 
performing the summation in Eq.([4.6|) (for generic k or d). To this end, we first apply the 
following formula to Eq.(^ 



KJz) = ^} Z / 2 } V r e - zt (t 2 - ly-^dt ■ Rev>-1, Rez>0, (4.7) 

and note that the k < 2 {d > 1) and k > 2 (d < 1) cases involve different calculations. Here we 
put a few remarks. For the convergence of the summation over n, we have to assume 

fj, < M (4.8) 
13 



in order to satisfy the condition 

\ e -Wt±P)\ = e -W*±ti < ! . t > ! . (4.9) 



The k = 1,0,-1 (d = 2,3,4) cases are calculated in Ref.[22], and the k = 3 (d = 0) case 
corresponds to Ref.[^], although they did not discuss this representation (Note that their 
argument belongs to our second representation). 



(i) In the first case, k < 2 [d > 1), the condition on v in the formula (O) is satisfied as it 
is; i.e., v = | — 1 > — ^, and the calculation is almost parallel to Ref.[^]. We then just write 
down the result 

where note that the minus sign of the denominator differs from the previous formula (Eq.(3.13) 
in Ref . |f22|| ) . For later convenience, let us derive another formula from this result. Performing 



the change of variable u = M(y/p 2 + 1 — 1), and applying the formula T(z)T(l — z) = 7r/sin7rz 
with z = {d — l)/2, we derive 

Of'(M) = -|r(|)M>-<= s ,n(^) f ^XS-'r * + <P ~ ~ P) ' (4 - H) 
where we have introduced the compact notation 



X±(p) = M^/p 2 + l±fL . (4.12) 

(ii) In the second case, 2 < k [d < 1), the calculation is similar to the case (i), excepting the 
point that we apply the formula (|4.7| ) for u = 1 — i instead of v = ^ — 1 owing to the relation 
K v {z) = K^ v (z). Then the summation over n in Eq.([4.6|) leads to the Lerch transcendental 
function: 

Of{M) = a|°° dt(t 2 - 1 )^ e -mt+P)^ e -mt+H) } d _ 2; i) + 

aT(3-d) f°° , , o /■(°+) (-z) d - 3 dz . , 

' dt(t 2 -l)— - { _>^ rA + (A -> -A) , (4.13) 



2vri 7i v Joo 1 - e *+/J(AfH-/i) 

where just for conciseness we have defined the coefficient 



r(^)V r(|)4 
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The z-integrand ( 4.13| ) possesses poles at z = and z = —(3Mt ± /3p, + 2niri (n G Z). Since 
—(3{Mt ± jti) is a negative value, there is no pole in the positive region on the real axis other 
than z = as long as /3u> + — is not an integer. With the change of variable t = \/p 2 + 1, and 
a replacement of one derivative d/dz by d/dp, we have 



1)1 dz k ~ 2 



fc-2 







V 



L(. 



i 



dp + (fi—* -ft) . (4.15) 



2=0 



dp V 1 - e*+/'- x '+ 

Because of fc > 2, the surface terms from a partial integral vanish in Eq. ( |4.15|) , and finally we 
find 

pk 3 



of\M) 



2T{ k±] L) (-l) 1 - k l3 k - 2 {k-2) 8 k ~ 2 



V 

1 _ e z+l3X + ' 



2 ' Af(Jfe-l)! dz k ~ 2 
More explicitly, we present the results for k = 2,3,4 as follows: 

1 



dp +(fi^> -p.) ■ (4.16) 



2 = 



Of{M) 



)(2) 
>(3) 



M 



0f\M) = £[ 



0f{M) 



1 _ e /3(M+A) 
(l_ e /3X + )2 

r 1 



dp + (/Z — ► - 



2M 3 



(4.17) 
(4.18) 
(4.19) 



. 1 _ eP(M+i£) (]_ _ e /3(M+/2))2 

For even values of fe, one may directly perform the summation ( f4.6|) without using the integral 
representation (4.7); e.g., using Ki/ 2 (z) = y / ir/2ze~ z and K s / 2 (z) = y/-ir/2z(l + z~ 1 )e~ z for 
k = 2 and 4. We will see in the next section that the results ( [4.11 ) and ( 4.17[ )-( ^.19[ ) can be 
reproduced from the second representation as well. 



5 The integral formulae from the second representation 



In this section, based on the second representation ( 0.14 ) case, we derive some more formulae 
on of(M), with re-deriving the results of the previous section. From Eqs.( |3.30| ), ( 0.36 ), and 



(4.2), we have the following form of the pure thermodynamic part to start with: 

W(M\= (4?r)1/2 I™ J^r-sW 2 -?) £ 

nez.^o 



Of (M) = M^ r<fa S ¥e-^-A 2 ) £ e-^) 2 cos(4nvr SA l//3) . (5.1) 



The summation can be converted to the integrals on the contour (see Fig. la): 



0f(M) = 1(4.)! f dss^e-^^- ^-y^dz + (£ -> -/2) . (5.2) 
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-•2 



After performing the s-integration, one may deform the contour C„ into another one C$ (Fig. lb) 



in the same way as Ref.[25], thus obtaining 



^ v ; v V2vri ioo e^ 2 "^) - 1 



<iz + (// — > — /i) 



(5.3) 



Parametrizing the circular part of C$ (centered at M) by z = M(l — 5e 1 ^), and using p 2 = z 2 — 1 
for the remaining parts of the contour, we arrive at 



p V 7 aAF V v 2 ; LA/25 - 1 



dp + 



,p{M-p) _ 1 



(5.4) 



Substituting k = 3 and e = 2, Eq.(IT3) reproduces the result (without background field) of 



Ref.|25|. 




Im z = 




C<5 



(6) 



Re/x = fj, 
(a) 

Figure 1: The contours and Cg. (a) The impounds an infinite number of poles 
at z = Jx + 2niri, where n denotes nonzero integers, (b) Coming from z = oo + ie and 
going to z = oo — ie, the C$ impounds the pole at z = M. 



Here are a few remarks on the computational difference between Ref.[25| and ours. We have 



extracted the pure thermodynamic part (5.1) from the beginning, and this fact is expressed by 
the circle at p, in C^. On the other hand in the method of Ref.[^], one needs an algebra to 
separate the non-thermodynamic part by using 

-1 1 



1 + 



(5.5) 
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The first term on the r.h.s. of Eq.(|5.5|) corresponds to the non-thermodynamic part, and the 
second term to the following expression instead of Eq. Q5.2p : 

ps(z— iui) 2 p p s(z—iu>) 2 



0g)(M) = (47r) 



dss 2 e 



-sM 



1 



1 



+ 



-f3(z—fi—iuj) 



-}dz , (5.6) 



2vri lJc + e P( z -P- ilZ >) 

where the contours C± run from /i db o — too to [i db o + too for a small value of o. Although one 
can of course obtain the same result Q5.4| ) from Eq. Q5.6| ), it is not obvious at first glance that 
the C- integral term corresponds to the /x — ► —p, term of Eq. (|5.3| ). 

Now, further results depend on the cases (i) and (ii). In the case (i), k < 2, the second term 
in the square brackets on the r.h.s. of Eq.(|5,4|) is proportional to 5 1 ~ k ^ 2 , and vanishes as 5 —* 0. 
Therefore Eq.( |5.4| ) reproduces the first representation result ( |4.11 ). 

In the case (ii), k > 2, the 5 l ~ k l 2 term diverges as 5 — ► 0. We expect that it can be 
canceled with the lower surface term from the first integral in the square brackets on the r.h.s. 



of Eq.( 5.4). This can be seen by a partial integral with applying the formula 

d l 

dp p p 2 ^p 2 + 1 

to the p- integration in Eq.(|5.4j) as follows: 



(5.7) 



oo pl~k^p2 _|_ n-1/2 



V2S 



dp 



^(2<5) 1 - fc /2 



fc-2 



(5.8) 



+ 



1 



V 



3-k 



(p 2 + I)" 3 / 2 M(3{p 2 + 



+ 



( e /3X + _ !)5 



}dp. 



The first term on the r.h.s. of Eq. (|5.8|) cancels the divergence as expected, and we derive 



Of\M) 



(5.9) 



P 



3-k f 



(p2 + l)"3/2 M ^ p 2 + !)-l e /3X + 



1 



+ 



( e /3X + _ 1)2 



}dp + (p -/i) 



Putting t = p + 1 in Eq.( |5.9|) , we finally obtain the following concise expression: 



0j fc) (M) 



4vrr(- M 



,_ fc sin(VT) 



(2-fc)7T A 



Vt 2 



—2-k d 

1 at 



t 



-1 



ol3(Mt+p,) 



dt+(fi^> -p) , (5.10) 



which is an alternative representation of Eq.( [4.16| ). As a result, we have derived two expressions 
in each case: Eqs.( 4.10 ) and ( 4.11 ) in the case (i), and Eqs.( 4~16| ) and ( 5.10 ) in the case (ii). In 
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Eq. fl4.16j ) we have to perform k — 2 derivatives (after one integration), while just one integration 
in Eq.( |5,10| ), whose integrand thus contains a common Boltzmann factor for all k (> 2). 

Let us check the consistency of Eq. ( |5.10| ) with the previous result flOp for k = 2,3, 4. The 
k = 2 case ( 4.17 ) is immediate from the representation (|5.10| ), however let us handle these three 
cases simultaneously. Using (p 2 + 1) _1 = 1 — p 2 (p 2 + 1) , we first divide the second term in 
the curly brackets in Eq.( |5.9| ), and then perform a partial integral on the r.h.s. in the following 
quantity: 

p d 1 



1 3 _ k Mpp 2 (p 2 + 1)-V X + 



3-fc 



V 



io (,e f '" + — i-)" Jo a/p 2 + 1 ®P eP x +-l 

Changing the integration variable by t 2 = p 2 + 1, we arrive at 



)dp 



4vrr(-)M 



-7T r 



(2 - k)rr 



where 



Vt 2 



1 



■2-fc rA; - 3 



N(t) - N'(t) 



dt 



and N'(t) = dtN(t) means the derivative of the function 

1 



N(t) 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



e f3(Mt-p.) _ i ■ 

For the k = 2 and 3 cases, this expression is sufficient to see the consistency in each. However, 
for the k = 4 case, notice that we have to extract a singularity, which cancels a zero from the 
sine function in Eq.( [5.ld ), and we hence perform a partial integral once more: 



1 



k 



4-k d 

M 



(k-3) 



N(t) 



t 



N'(t) 



dt 



2<k<A 



(5.15) 



where the surface terms can be dropped only when k < 4 (also in Eq.( |5.1ll )), and this is the 
reason why this formula is valid for 2 < k < 4 (—1 < d < 1). Substituting k = 2,3, 4 — e 
(with e — > 0), we verify that these expressions for reproduce Eqs.( [Q7l) -( [4.19| ), In this sense, 
Eq.fljuj]) is equivalent to Eq.(gl|). 

For further values of k, one should repeat the similar calculation for each interval between 
zero points of the sine function. For example, for k = 6 and 8, we obtain 



of(M) 
of{M) 



7T 



16M 5 



8M 7 



3iV(l) - 3iV'(l) + N"(l) 
15JV(1) - 15iV'(l) + 6N"(l) - JV"'(1) 



(5.16) 
(5.17) 
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These pure thermodynamic 'partial' amplitudes exactly correspond to the N = 4 and 5 pure 
thermodynamic amplitudes in D = 3 when Vqu = V (the first category). In D = 4, odd k 
integers are the similar cases. For the second category; for example, the iV-photon amplitude 
in D = 4, we need to combine the quantities from k = — 1 to 2N — 5. These cases can not get 
rid of non- integrated quantities such as Eq.( [4.18 ). After all, Eqs.(4.16) and ( 5.1 0| ) describe the 



general amplitude formulae which contain all Feynman diagrams. 

6 The (3 -> oo limit 

In the above arguments, we have assumed /x < M, and all [5 — ► oo limits vanish because of 



it: 



Of ^ 0. (6.1) 



To obtain a nontrivial (nonzero) limit, we should remove this condition after all. In this sense, 
the //-dependence of the master amplitudes is non-analytic. To see this, we need to transform 
the function of in a form indicating a Bose (or a Fermi) distribution, and we already derived 
this kind of representations in Sections || and |5[ 

Let us begin with the case (i) k < 2; in particular the case of k = (which is also an 
odd dimensional case D = 21 + 1 by the way). From Eq.( [4.10|) with changing the variable 
E = u + M, we have 

0f{M) = -2^ £ l± %^ + (M' - V) , (6.2) 
where the plus/minus signs correspond to the Fermi/Bose statistics, and we put 

fi' = (j, — iui . (6-3) 



Eq.( |6.2| ) can be interpreted an effective action or total energy density with the chemical potential 



fj,' and the mass M ( < //, E ). Now, we understand in Eq^^) that the becomes zero as 



P — > oo if n' < E, while the takes a nonzero value if E < //. The similar arguments apply 
to the case of generic k value in the following way. Since we observe 

e m-P+u) = e f3(u+M-^ e iP(o+^) P^c Q { OI u<fl _ M , (6.4) 
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we have the finite upper boundary on the integral ( 4.10| ) at u = /j — M, thus obtaining 



Og\M) = Jim Of (M) = r M (u 2 + 2Mur k ' 2 duO( l JL - M) (6.5) 

P^oo ' T( £ 2 ± ) JO 

~ 2 V^ f r , 2 ~ k 2-k 4-k M - /I 



r (4^fc)V ; ; v 2' 2 ' 2 ' 2M ; ^ ; ' 

The explicit results for fc = — 1, 0, 1 are 

0«(M) = -2axccosh(/x/M)0(/z-M) , (6.6) 
0®(M) = -2V5F(/i-M)0(/i-M), (6.7) 



O^^(M) = -2 / u v / / u 2 -M 2 + 2M 2 arccosh(/i/M) 6(fi-M) . (6.8) 

In one of the = cases (D = 3, N = 0, I = 1, e = 2) ||, Eq.Q coincides with an exact 
result in ]27]| . 

The case (ii) k > 2 can also be estimated in the same way. Applying the similar treatment 
as above to the Boltzmann factor in Eq.( 5,lC| ), we then obtain the nonzero limit given by 



0«(M) = -v^)M^!g^ . (6.9) 

Performing the i-integration, some of explicit results (the k = 2, 3, 4 cases) can be shown as 
follows: 

Og){M) = ^^-l)e(n-M), (6.10) 

OgKM) = -^\j^) 2 - l6 ^- M ) ^ ( 6 - U ) 
0«(M) = --^Le(^-M). (6.12) 
The calculations are rather straightforward for further k values, and the general formula (|6 

(k) 

suffices. The e dependence is gone away from the Ooo ■ This is natural because there is no 
difference in the kinematical factors (hence in 'partial' amplitudes) between boson and fermion 
loops at zero temperature. 



7 Conclusions 

In this paper, we studied the thermodynamic generalization of the one-loop amplitudes 
which can be cast into the BK master formula ( |2.9[) at zero temperature (with zero chemical 
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potential). We followed the two-step procedure: first, calculating the path integral S 1 sum- 
mation ( p.l| ) to introduce the temperature, and then performing the shift manipulation ( |1.3| ) 
to insert the chemical potential in a loop. This procedure has been applied parts by parts to 
the kinematical factor /C, the normalization factor J\f and the vertex structure function (re- 
duced kinematical factor) V, and we have derived the general formulae for these parts. From 
Eqs. fl3.31 ), ( 3.34| ) and fl3.35|) , the thermodynamic Appoint amplitude of general form is thus 
summarized as 

T N=c E a in^ / da x / da 2 --- / da N Af" . (7.1) 



One can check the validity of this master formula in various cases; the free effective potentials 
and the photon polarization in Appendix B, the effective potential in a constant magnetic field 
in Ref.(|§, and the vr° -> 2 7 decay in Ref.||] etc. 

The detail analyses on the building blocks Vp^ and KLp^ have given us useful information. We 
have realized that the total thermodynamic kinematical factor ICp^ behaves as the thermody- 
namically generalized normalization to the zero-temperature type correlator shown in Eq.( |3.17j ). 
This is certainly a useful result and makes calculations for large N much simpler than ever; 
namely we have only to attach the new normalization A//?^ to the Appoint scalar correlator (with 
switching — ► o^.). Another interesting point is that the pure thermodynamic part of a vertex 
structure function can be expressed in terms of either vanishing or ^-differential operators, and 
this fact makes us possible to decouple the master formula ( |3.31| ) into T^ 1 = Tn + fjy in a 
nontrivial way (cf. Eqs.( |3.31 ) and ( 3.33| )). It is worth noticing that we do not decouple the 
vertex structure function but the kinematical factor only. From Eqs. (|3.33| ), ( 3.34| ) and ( |3.36| ), 
we then conclude that the pure thermodynamic master amplitude is simply given by 

Apart from r^r, the pure thermodynamic part is renormalization free, and hence Eq.( |7.2| ) is 
the final form to apply our integral formulae derived in Sections |J-||. More explicitly, we have 
(from Eqs.(j3~T3D and (gg)) 



^ = c(4vr)^ Y, "Ini-iT-^ f^l f da 2 ■ ■ ■ f da N of +l ~ D) {M) a . (7.3) 

~!„ OUT Jq Jo Jo H m J -*m 2 +6 ! 
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and in particular for the [3 — > oo limit with /j ^ 



lim f^ = C (4vr)^^a /0 da x f da 2 ■ ■ ■ [ da N Og l+1 ~ D \M) 

'-♦oo leZ J JO JO 



m' 2 —*m' 2 +bi ^ 



The thermodynamic master amplitudes (7.3) and (7.4) are provided with Eqs.( 4~l~6| ), ( |5.1C| ) and 



(k) 

and several explicit 'partial amplitudes' O ^ are also given up to k = 8. 
In the master formulae ( ff^ ) and (|7.3|) , the 'partial' amplitudes or 0„ (s-integrals of 
the pure thermodynamic kinematical factor) are the fundamental computational blocks in our 
formalism, and in Sections || and [|, we focused on some mathematical aspects and techniques 
how to calculate these integrals. We derived various representations and formulae on these 
integrals for arbitrary values of N, D and e. It is obvious that a variety of equivalent formulae 
makes easy to find some relation and consistency with the results obtained by different methods. 
These should generally be related by analytic continuation, and a different analytic expression 
is certainly useful by definition; i.e., another formula can cover the region which can not be 
reached in an original representation, and it also overcomes some defect of an inconvenient 



representation. For example in Ref. [22|, a dual rotation is used to obtain an electric gap 
equation from magnetic one in a four-fermion model. We expect such kind of utility when our 
formulae are applied in more explicit stages. 

Although we could have possibly simplified our results furthermore through a certain tech- 



nique [23], our formulas were sufficiently convenient to extract the nonzero values of zero tem- 
perature limit with /i kept finite. For this purpose, it is necessary to have the integral represen- 
tation, such as Eqs.( 4.10 ) and ( 5.10[) , which clearly indicates a non-analytic cut in its integrand 



at (3 = oo. We also had to go beyond the condition [i < M in these representations, however 
this might be justified by analytic continuation, and at least we know that the k = case 
coincides with the exact result. 

We have examined the model independent structure of the thermodynamic BK master 
formula with consulting several simple examples. The model dependence appears in the vertex 
structure functions, and hence one has to evaluate the vertex structure in the first place for 
explicit calculations in each model. This task is case by case and will become more lengthy as 
TV increasing for the cases involving more bosonic x-fields in the vertex operator Vj parts such 
as photon, gluon and pseudo-scalar particle. However, our systematic prescription is certainly 
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promising to obtain the (pure) thermodynamic iV-point amplitudes in a straightforward way as 
long as the belongs to the general form (3.34). Finally, we should not forget the advantage 
that the world-line formula encapsulates all necessary Feynman diagrams in a single integrand. 



Appendix A. Insertion of chemical potential 

In this appendix, we give a brief explanation of the shift procedure (|3.9|) . First consider 
the vacuum case Q = 0. In this case, M 2 = m 2 and 

- Mn(* + m 2 ) = \ £ / 7^1 f . (A.l) 



P^J (2tt) d - 1 J s 
We usually perform the insertion of chemical potential in terms of the shift 

P 2 = Po + f -> K + in? + f (A.2) 
with the internal Matsubara frequencies 

27T 1 , A N 

w n = — (n+-), (A.3) 

and we have 

a3^ /- 2fl A) • 47 i" 



= e-«e 3 (^,M^) , (a.4) 



where we have defined 



27F * I A r\ 

Mo = A» - —g • ( A - 5 ) 



Thus, we prove Eq.(3.14) for A r = 0: 



Eq .(A.l) = i jT £ ( 4„)^« ! *«- > )e,(^, . (A.6) 

For further nonzero values of N, the proofs are straightforward, and we shall not explain the 
details anymore. For example, one can find the N = 2 case in Appendix B. 
Instead, we add a supplemental interpretation. Applying the transformation 

^^e 3 (^,4) = es(^,A , (A.7) 



P K P ' I3 2 ' Ay 2vr ' 4vrs ; ' 

we rewrite Eq. (|A.6|) as 

Eq .(A.6) = X" = f -(<M^-"" 2 e 3 (if^,#) . (A. 

JO S Z7T 47TS 
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Here p,Q now appears only in the first argument of the G3, and remember that ui also appears 
in the same place as fio does for N 7^ 0. Taking account of exponent's additivity in the first 
argument of 03 ( q. v. Eq. (|3.5| ) ) , one can imagine that the final result is given by the replacement 
of p,Q with fiQ — iLJ ( = /2 ). This result coincides with the /Cg„ obtained by the shift manipulation 
in JCp. 



Appendix B. Photon self-energy part from Feynman rule 

In this appendix, we rearrange the N = 2 Feynman amplitude of photon scattering into the 
world-line formula at finite (5 and [i. The photon self-energy part (in D = 4) by the Feynman 
diagram technique is written in the form 

where p^ = (u n ,p), kff = (uj^, and = e* i = 1, 2. We also make use of k^ = k^ = —k^. 
The QED coupling e is set to be the unity for simplicity. For convenience, we decompose the 
11^ into the following three parts: 



in terms of 



where 



Tr[A (/ - MW] = W 1 + W 2 + W 3 , (B.3) 



W X = 2ei-e 2 (p 2 + {p - kf) , (B.4) 

W 2 = 2(ei • e 2 k x -k 2 -ex- k 2 e 2 ■ h) , (B.5) 

W 3 = -2ex ■ (2p - h)e 2 ■ (2p - h) . (B.6) 

Applying the Feynman integral formula 



1 



dss I 1 d ue skHu-u*) e -s(p-(l-u)kf j ( B 7 ) 



p 2 {p-k) 2 Jo Jo 
we easily rewrite IT and II2 in the following forms: 

Hi = ^ E r 77^72 Aie-*"^ 2 )e-^+«) 2 - „) , (B.8) 
P „ Jo {4irsr' z Jo s 
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H 2 = ( Cl • e 2kl ■k 2 -e l . k 2 e 2 • ky)l £ jf jf * ^-^(^e-KW , (B . 9) 

where 

a = (1 - u)o; fc • (B.10) 

The II3 can be rewritten in the similar way by using Eq.( |B.7D and shifting p — ► p + (1 — u)fc in 
the p-integral. After some algebra, we have 

- ^r^->-»%--<^ (B.11) 

[ -(ejeg - ei • e 2 ) - { 2ejK + a) - (1 - 2u)ei • fc }{ 2eg(w n + a) - (1 - 2u)e 2 • k } 



x 

s 



Gathering Eqs.flBj), Q and (pUD , we obtain 



P P^Jo (4^)3/2 io 



n 

r 2 

x ~ei • £2! 5(1 — u) — 1 } + (ei • e 2 &i • /c 2 - ei • /c 2 e 2 • k\ ' 
L S 

- (1 - 2n) 2 ei ■ ke 2 ■ k — 4e^e 2 l {(w n + a) 2 - — } 

+ 2(l-2u)(eJe 2 -A; + e 2 l e 1 - fe)(a? n + o)] . (B.12) 

Now, we can eliminate the forms (u> n + a) m ; m = 1,2, in the summand of Eq.( |B.12^ ) due to 
the following operations: 

" T 8a £ e ~ S( " n+a)2 = E K + a)e- s ^ +a ) 2 , (B.13) 

^ n n 

(^<9a) 2 £ e ~ sK+a)2 = E[K + a ) 2 -^] e " S( " n+a)2 • (B.14) 

The a-dependence of the summand exponential can be transformed into a linear exponent form 
by using the Jacobi transformation: 

^ ^(_i)« e -4f i e -^ . (B.15) 



V47TS 

Then we are allowed to perform the following replacement in the transformed summands 

^ -> (B.16) 
25 



and we finally derive the formula 



U = 2 f°° , d Z,o C due- sk2{ - u - u ^Y{-l) n e-^e in ? a 
P Jo (4vrs) 4 /2i ^ > 

r i 9 n 2 /? 2 in 6 . , . , 9 

x ejeg— £- + — (2u-l)(eSe 2 -fci-egei-fc2) 

2 

+ -ei • e 2 { <5(1 — -u) — 1 } + (ei • e 2 /ci ■ fc 2 — ei ■ fc 2 e 2 • k\) 
s 

+ (1 - 2ufei ■ k 2 e 2 ■ h ] . (B.17) 

For /i / 0, we have only to make the standard shift ([O]), and this causes the following 
change of the a defined in Eq.(B.lO): 



a -> (1 - u)u k + ifj, . (B.18) 

This modification exactly corresponds to the shift (|3.S|), Note that the a coincides with u> with 
fixing u\ = 1. 

Appendix C. World-line method for photon self-energy part 

In this appendix, we illustrate how we obtain the thermodynamic version of the vertex 
structure function in the case of N = 2 photon scattering. In the world-line formalism, the 
photon self-energy part at zero temperature can be obtained from the formula 

T 2 = --2T / -e~ sm & VxVij exp[- / {-Mr) + -^drjV^ , (C.l) 
2 Jo s J Jo 4 2 

where Vj, j = 1, 2, are the photon vertex operators 

Vj = -ie [ S drjiej ■ x + 2# ■ erf ■ kj)^)^^^ , (C.2) 



o 

and x = d T x etc. Here we follow the standard world-line notation r, which is related to the 
main text notation a by 

Tj = soj . (C3) 

Eq.( C.l ) is known to become 

/rs rs 
dn 
o Jo 



1 d r°° r/s 

r 2 = --2T / — / dn I dr 2 VxIC (C.4) 



'We set the QED coupling e = 1 in the following (as in Appendix B). 
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with the vertex structure function 

V = ex ■ e 2 G% + £l • k 2 e 2 ■ ki{G l g) 2 + (ei • e 2 h ■ k 2 - e x ■ k 2 e 2 ■ k^G^f , (C.5) 
and the kinematical factor 

= (4 7rs )-^/2 e -sm 2 e fc 1 .fc 2 Gi ; 2 ; ( C _g) 

where 

G 1 i = \r 1 -T 2 \- {7l ~ T2) \ G¥ = sign{ n -T 2 ) . (C.7) 
s 

The path integral normalizations are chosen to be 

j> Vxe~ /o = (4vrs)- D / 2 , j> V^e' f° ^ dr = 1 . (C.8) 

It can be said that the kinematical factor is defined by the insertions of the 4? scalar vertex 
operators 

Vj = f dTjexp[ikj ■ x(r,)] ; j = 1,2 . (C.9) 
Jo 

According to the program presented in Section ||, we are led to calculate the following path 
integral at finite temperature: 

Id ,. so r°° ds _2 r , r r s ,i .„ 1 



= -_2T(ie) 2 / _ e - sm2 i exp[- / (^i 2 + -V • 

2 Jo s J Jo 4 2 

x (II [ dTie^^Y^i-^Tz'^ z in ^ {Tl ~ T2)lJJk 



x (e^ + e 1 •i; + 2^-eiV'-fci)(n)(eo— + e 2 • i + 2# • £2^ • ^)(-r 2 ) . (CIO) 
s s 

Using the Wick contraction method with the correlators 

< ^(nK(T 2 ) >= -<TGi, 2 , < ^(n)^(r 2 ) >= \sTGf , (C.n) 

one can verify the coincidence of the with the Tig derived in Appendix B; we have arrived 



at the form 



n, = I 2 f T-e— 2 (f[ /%)£(-irX^(V + V) , (C.12) 
2 Jo a Jo 
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where 



V = e e — r~ H ( e o £ 2 • fci - e £ i • k 2)G 



B i 



(C.13) 



and the is the nth mode of the bosonic two point function at finite (3 defined by 



x(Tj) 



i=i 

n 2 8 2 - a , \ i 



(47TS) 



-D/2 e ~k 2 G\ 



(C.14) 



Here Eq.(C.12) with ( |C. 13 ) and (C.14) reproduces the Feynman rule result ( B.17| ) with rescaling 
Tj = scij and fixing o\ = 1 with o"2 = u. Note that G^ behaves as 2u — 1 in this respect. 

Let us further rewrite the above V in an operator form suitable to Eq.( |3.3i ). In the following, 
we shall not employ the fixing a% = 1. As shown in Appendix B, we can replace in(3 with d/du), 
where 

U}=((T l - a 2 )uj k , 



and we therefore find 



1 /"OO ^ P 1 

rf = -2T / d ss (n / ^)V s x/C 

2 Jo jLj Jo 



with having 



and 



= E(- 1 ) n 4 n)e 



(C.15) 
(C.16) 
(C.17) 
(C.18) 



where the change of variables = soi in G^ should be understood. Since we already justified 
the shift (|3.9|) in the end of Appendix B, we can use the following parts for the finite \jl case 



d d 



(C.19) 



Substituting these in Eq.( |C.16 ), we obtain the two-point function Ig > which therefore coincides 
with the corresponding Feynman rule result suggested in the end of Appendix B. 
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